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In this work we study the effect of thermal processing of exfoUated graphene on mica with respect to changes 
in graphene morphology and surface potential. Mild annealing to temperatures of about 200 °C leads to the 
removal of small amounts of intercalated water at graphene edges. By heating to 600 °C the areas without 
intercalated water are substantially increased enabling a quantification of the charge transfer properties of 
the water layer by locally resolved Kelvin probe force microscopy data. A complete removal on a global scale 
cannot be achieved because mica begins to decompose at temperatures above 600 °C. By correlating Kelvin 
probe force microscopy and Raman spectroscopy maps we find a transition from p-type to u-type doping of 
graphene during thermal processing which is driven by the dehydration of the mica substrate and an 
accumulation of defects in the graphene sheet. 

Because of its exceptional properties, graphene, an atomically thin layer of carbon atoms in sp^ configuration, 
has attracted great interest in the scientific community. Especially its extremely high thermal and electronic 
conductivities, mechanical strength and impermeability to even the smallest gaseous molecules make 
graphene a highly interesting material for applications in e.g. electronics, photovoltaics, sensors and water 
desalination to name only a few^"^. 

Being a material which is nothing but surface, graphene is known to be largely influenced by the 
underlying substrate and the environment. Pure graphene should have no intrinsic charge carriers but in 
practice it often shows p-type doping caused by a combination of air humidity and oxygen molecules^"^^. 
Excess charge carriers (electrons as well as holes^^"^^) can also be transferred from the substrate, which thus 
plays a major role for the performance of graphene devices. In addition, charge puddles and impurities in 
the substrate may act as efficient scattering sites, as observed for graphene on Si02^^'^^. It is therefore of 
utmost importance to choose the substrate carefully and to quantify its specific doping and charge transfer 
properties. One of the most promising substrate candidates for future graphene devices is mica^^ Using 
molecular beam epitaxy it has been shown that it is possible to directly grow graphene on mica^^, which is an 
important factor for industrial applications. Mica is a layered silicate material that exhibits a very flat surface 
after cleavage. As a consequence, exfoliated graphene is also extremely flat with a root mean square (rms) of 
d = 24.1 pm^° measured by atomic force micoscopy (AFM). However, this flatness can only be achieved by 
exfoliating graphene in a dry atmosphere. Exfoliation in ambient leads to intercalation of water layers^ ^"^^ 
which could prove detrimental for electronic devices. Raman spectroscopy, Kelvin probe force microscopy 
(KPFM) and scanning tunneling spectroscopy performed in ambient revealed a p-type doping from the 
substrate, that is effectively blocked by the intercalated water layers. It was suggested in ref. 24 to actually use 
water to modulate the electronic properties of graphene, but in practice controlling water islands locally 
would be extremely difficult to realize. Therefore one is confronted with the problem to remove the 
intercalated water. 

In this work, high temperature processing in ultra high vacuum (UHV) of graphene on muscovite mica 
has been performed in order to study the effect of heating on the interfacial water layer. To this end, 
pristine and thermally processed samples have been thoroughly investigated with respect to morphology 
and work function. We present data from in situ non-contact AFM, showing that basically a complete 
removal of water cannot be achieved because decomposition of the mica substrate limits the maximum 
temperature. We therefore proceed by recording changes of the doping level of graphene due to in situ 
heating with locally resolved KPFM to probe the effect of intercalated water on the charge transfer 
between graphene and its substrate. This data is complemented by Raman spectroscopy measurements, 
indicating a clear correlation between thermally induced structural changes and changes in charge carrier 
concentration. 
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Figure 1 | (a) Non-contact AFM topography of a pristine graphene flake on mica. The SLG flake is decorated with water islands, (b) Zoom-in 

on the white box marked in (a). The red line profile denotes the height of the SLG flake on mica while the black line profile denotes the height of a water 

island, (c) Surface potential map corresponding the the topography in (b). Water islands lead to a decrease in surface potential by up to 200 mV. 



Results 

For our studies we used mechanically exfoliated graphene flakes 
supported on muscovite mica. In fig. 1(a) a non-contact AFM topo- 
graphy of a typical graphene fiake on mica after introduction to UHV 
conditions (phase ^ 1 X 10~^° mbar) is shown. The zoom-in topo- 
graphy image in (b) shows that the single layer graphene (SLG) fiake 
is decorated with small islands. The red line profile in fig. 1(b) 
denotes the height of the SLG on mica which is about 0.4 nm. 
Whfle this value corresponds wefl to the nominal interlayer spacing 
of graphite, we observe that the measured heights frequently vary 
between 0.4 nm and 0.9 nm. This is in agreement with values re- 
ported for SLG on mica varying between —0.8 nm and 1.0 nm^^'^^. 
By exfoliating graphene on mica in ambient conditions, water layers 
wiU be trapped at the interface between graphene and mica which 
results in increased step heights measured by AFM^^'^^. These inter- 
facial water layers, IWL cafled in the foflowing, have been found in 
different thicknesses from one to three IWL. The smafl islands which 
are up to 1.2 nm in height (see black line profile in fig. 1(b)) are 
attributed to water on top of graphene, as they can be easily removed 
by mfld heating to 180°C, which is significantly below the temper- 
ature needed to clean graphene from contaminations^^'^^. They can 
easily be identified in the surface potential map as well (see figure 
fig. 1(c)). Here, the water islands decrease the surface potential of 
graphene by up to 200 meV. 

To obtain quantitative data from our KPFM measurements, the 
work function of the AFM tip is calibrated on few layer graphene 
areas on the sample. The work function value of graphite can be 
assigned to few layer graphene areas where the surface potential 
depends on the layer number^^. This procedure yields ^flg = 
^Graphite ^ 4.65 eV for few layer graphene (FLG). With the tip cali- 
brated in this way, the work function of single- and bflayer graphene 
can be calculated using the relation — ^Graphite + ^(CPDpiQ — 
CPD^g)^ whereas CPD is the contact potential difference measured 
by KPFM between the AFM tip and the sample. In graphene, a 
variation of the work function is directly correlated to a shift of the 
Fermi energy /S.Ep with respect to the Dirac point^°. Based on the fact 
that the work function of graphene is around 4.57 eV^^'^^, a work 
function higher than this value can be attributed to />-type doping 
and vice versa n-type doping for a lower work function. 

Next we tried to apply a thermal treatment to remove the inter- 
calated water film. To this end, several heating steps where per- 
formed at different temperatures and exposure times. The pristine 
SLG fiake in fig. 2(a) with the continuous IWL exhibits a work 



function of 4.83 eV which is lowered to 4.66 eV in areas where water 
islands are adsorbed on top of it as seen in fig. 1(c). Heating this 
sample to about 180°C for 1 hour can be used to remove the adsorbed 
water islands but not the interfacial water as can be seen in fig. 2(b). 
As the intercalated water cannot penetrate the graphene sheet, only 
small areas at the edge of graphene are freed from water. Because 
grain boundaries and defect sites are very rare in exfoliated graphene 
sheets, no effective alternative desorption pathways exist. The work 
function map shows a simflar work function distribution, the work 
function is decreasing with increasing layer number. However, the 
work function of graphene in direct contact with mica is about 
0.17 eV higher compared to graphene lying on the interfacial water 
layer. This is an important finding as it suggests that pristine gra- 
phene and mfldly heated graphene both show p-type doping. In 
contrast to KPFM measurements in ambient (where p-type doping 
is always present because of humidity), this can be taken as proof that 
graphene on mica is indeed p-type doped by a charge transfer from 
the mica substrate. 

The last images in fig. 2 show results that were obtained after 
heating the sample to about 600°C for 19 hrs. The area where the 
intercalated water layer has been removed has increased (dark blue 
areas in fig. 2(c)) but in large areas of the graphene fiake intercalated 
water can stifl be clearly seen (light blue areas in fig. 2(c)). The work 
function has changed significantly over the whole fiake: From 4.8 eV 
(greenish area in fig. 2(b)) to 3.93 eV (bluish area in fig. 2(c)) where 
the water has been removed, and from 4.65 eV to 4.28 eV in the areas 
where the water could not be removed. In addition, new character- 
istic features have appeared (marked with "nanoblisters" in fig. 2) 
with typical diameters of about 10 nm and heights of 0.3 to 1.5 nm. 
Note that these features are quite different from the water islands 
found before heating (fig. 1), because they are much higher and are 
located underneath the graphene and not on top of the graphene 
sheet. We repeated the experiment with other fiakes and found that 
the observed changes are indeed characteristic but with small devia- 
tions with respect to absolute values (see fig. 4 for example). 

At this point we want to stress that the observed change in work 
function upon excessive heating is not only surprisingly pronounced, 
it also represents a transition fromp- to n-type doped graphene. Such 
a transition has never been observed for graphene on mica. Of 
course, any KPFM data has in general to be treated with great care 
due to uncertainties with respect to calibration. Therefore, to ensure 
that the work function values measured here are indeed correct, we 
undertook extra efforts to exclude the possibility that our calibration 
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Figure 2 | (a) Topography and work function map of pristine graphene on mica. The work function decreases with incresasing layer thickness. The sketch 
illustrates that a continuous IWL is beneath graphene while adsorbed water islands are on top of it. (b) Topography and work function map after heating 
the sample to 180°C for 1 hour. The water islands on graphene are removed and the IWL is not completely continuous anymore and graphene is partly in 
direct contact with the substrate. The work function is still decreasing with increasing layer thickness, (c) Topography and work function map after 
heating the sample to 600°C for 19 hours. The size of the IWL is further decreased and nanoblisters beneath the SLG flake are formed. The contrast in the 
work function map is inversed, as the work function is now increasing with layer thickness (Scale bars in all images are 400 nm.) 



material, few layer graphene, is affected by residual water (see suppl. 
material). Having excluded artifacts due to calibration, we must 
attribute the observed transition to heating-induced changes. 



To test whether the original state of the system with respect to 
work function can be restored, we exposed a sample to ambient 
conditions for three days after it had undergone the same treatment 



Pristine 



5.1 eV Heated/Defective 4.7 eV Ambient 



4.9 eV 




4.3 eV 



3.2 eV 



4.4 eV 



5.4 

5.2 F 

5.0 - 

i 4.8 P 

o 4.6 
c 



4.4 - 



I 4.2 I- 
4.0 - 
3.8 - 



I 



■ SLG 

• SLG+IWL 

A BLG 

▼ BLG+IWL 

I 



Pristine 



Defective 



Ambient Exposure 



Figure 3 | (a) Work function map of a pristine mica supported SLG, BLG and FLG flake. The work function decreases with increasing layer thickness, 
(b) After heating to 600°C the work function is inversed and intercalated water causes an additional work function contrast on SLG and BLG. (c) Work 
function map after three days of exposure to ambient conditions. The work function is again inversed, restoring the work function similar to its pristine 
condition, (d) Summary of the work function measurements from (a-c) (Scale bars in all images are 400 nm.) 
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Figure 4 | (a) AFM topography of a graphene flake after in situ heat treatment. Nanoblisters and areas with flat plateaus (ruptured nanoblisters) are 
marked in the image, (b) Raman spectra of pristine and heat treated single layer graphene on mica, (c) Raman mapping of the D/G peak intensity ratio of 
the graphene flake shown in (a). 



as described above. The results of this experiment are summarized in 
fig. 3. The work function of SLG in its pristine condition 5.09 eV (p- 
type) decreases to 3.90 eV (n-type) after in situ heat processing. The 
same transition can be observed for bilayer graphene (BLG) where 
the work function is decreasing from 4.88 eV to 4.35 eV. After the in 
situ heat processing, the sample was exposed to ambient conditions 
for three days. The work function map in fig. 3(c) shows a KPFM 
measurement after this exposure to ambient conditions. Again the 
SLG and BLG sheets exhibit p-type doping with a work function of 
4.82 eV and 4.78 eV, respectively. Thus, the initial values are almost 
completely restored. We attribute this to the exposure of water and 
oxygen present in ambient conditions, which causes a large p-type 
doping due to a redox reaction 

To find the origin of the observed contrast inversion in the surface 
potential, we investigated a sample in ambient conditions using AFM 
and Raman spectroscopy. In fig. 4(a), an ambient AFM topography 
of a similar treated - heating to 600°C in UHV - mica supported 
graphene flake is shown. The few layer part of the graphene flake is 
covered with nanoblisters that are up to 15 nm in height and about 
400 nm in diameter. Additionally, round plateau-like structures of 
similar diameter with a constant height of about 1 nm are observed. 
These flat plateaus could very well be the remnants of ruptured 
nanoblisters. Shown in fig. 4(b) are the Raman spectra of SLG gra- 
phene before (black line) and after (red line) heat processing. In its 
pristine condition the Raman spectrum consists of the well known 
first order non resonant G mode at 1588 cm"^ and the resonant 2D 
mode at 2645 cm~\ With a full width half maximum (FWHM) of 
24 cm"^ for the 2D mode, this Raman spectrum can be unambi- 
giously assigned to SLG^^ After heating in UHV two new modes 
have evolved in the Raman spectrum, the D and D' mode, which 
are both induced by disorder^^"^^; in addition the 2D mode is sup- 
pressed. The presence of the D and D' mode indicate a very high 
defect density. Raman spectroscopy can be used to quantify the 
amount and even the nature of disorder in graphene as well by 
determining the IdUg and IdUd' intensity ratios. Here, we find Id/ 
Ig= 1.7 and Id/Id' — 6.3, respectively. From comparison with pub- 
Ushed data we estimate that the defects here are mainly vacancy 
defects^^ and the mean distance between defect sites is on the order 



of Lo < 5 nm^^'^^. In an attempt to correlate the measured morpho- 
logy with the defect distribution we used Raman mapping. In fig. 4(c) 
the D/G-peak intensity ratio distribution of the measured sample 
area can be seen. Due to the very different lateral resolution of the 
two methods, a one-to-one correlation is not possible. Nevertheless, 
at many of the locations where large nanoblisters and plateaus (most 
likely ruptured nanoblisters) where observed in the AFM image, the 
Raman map does show areas with an increased D/G ratio (marked in 
fig. 4(c)). This strongly suggests that the formation of nanoblisters is 
accompagnied by a significant defect formation in graphene. 

Discussion 

The experimental evidence presented above suggests that the work 
function of graphene on mica is governed by different mechanisms 
which can be controlled by thermal processing. In the following we 
will discuss three possible mechanisms induced by heat: removal of 
adsorbates, substrate changes and defect formation. 

Any kind of adsorbate is expected to have a large influence on the 
properties of graphene. For example, graphene field effect transis- 
tors of graphene exfoliated on Si02 usually show very poor trans- 
port characteristics without proper cleaning by e.g. current 
annealing or annealing in forming gas of the device prior to mea- 
surement. Here, we are dealing with different adsorbates, namely 
water (from the humidity always present in ambient conditions) 
and air (O2, N2, etc.). Their effect on (non-defective) SLG has been 
investigated by several groups^'^°'^°"^^. Hitherto, only p-type doping 
of graphene on mica has been reported, but all data has been 
obtained under ambient conditions'^, where electrophilic adsorbates 
can be expected to mask the intrinsic properties of graphene. In 
general, the adsorbates obviously result in an effective p- type doping 
as confirmed by our data. As we could show, it is basically imposs- 
ible to completely remove the water from in between graphene and 
mica, but heating to 600 °C ensures that no water is adsorbed on top 
of graphene and that the intercalated water layer has been succes- 
fully removed at least in some areas. From our locally resolved 
KPFM data, we can easily resolve the difference between areas with 
the IWL present and areas without the IWL. This enables us to 
reveal the different and seemingly adverse effects that water may 
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have on work function measurements: With samples prepared 
under ambient conditions, intercalated water leads to a decrease 
of the p-type doping level (see fig. 3(c)), while for heated graphene 
on mica, intercalated water decreases the number of excess elec- 
trons, i.e. lowers the n-type doping level (see fig. 3(b)). These find- 
ings clearly underline the importance of UHV measurement when 
performing quantitative KPFM studies to determine intrinsic prop- 
erties of graphene. 

The second mechanism that may give rise to changes of the work 
function of graphene is the charge transfer from the substrate. It has 
been shown that the doping level of graphene is very sensitive to the 
mica surface configuration^^. From our KPFM data we find that the 
mica substrate does indeed show a change of its work function from 
4.6 to 3.4 eV upon heating (see fig. 2), indicating that the mica 
surface undergoes strong changes towards an electro negative sur- 
face. This transition is most likely driven by a depletion of potassium 
ions on the mica surface^l The fact that we observe a continuous 
change in the graphene work function upon heating from 5.1 eV (p- 
type doping) to 3.9 eV (n-type doping), can thus be interpreted in 
terms of a an increasing electron charge transfer from the mica 
substrate to graphene due to the increasing dehydration of the mica 
surface. 

Apart from direct charge transfer from the mica substrate itself, 
the heat induced changes of the mica susbtrate may have an impact 
on the work function as well. Significant defect formation in gra- 
phene can be detected after heating to about 600 °C where the 
dehydration of mica sets in^^'^^ Based on the observations presented 
above we therefore suggest that the following mechanism is giving 
rise to defect formation in graphene: Upon excessive heating mica is 
dehydrated because the chemically bound water evaporates from 
within the silicate. This water cannot permeate through the gra- 
phene layer and will thus accumulate and form nanoblisters trap- 
ping the water molecules. This is corroborated by the fact that 
graphene nanoblisters are reportedly accompagnied by defect 
formation^^"^^. In addition, in the process of the dehydration of 
mica, O2 is formed and thus oxidative etching as observed by 
Yamamoto et al. of graphene on mica may play a role for the defect 
creation formation as well^^'^°. 

In summary, we have shown that for mica supported graphene 
heating in UHV to temperatures of 600°C is no feasible way to free 
graphene of interfacial water layers which are confined during sam- 
ple preparation in ambient. The thermal processing of mica in UHV 
causes the formation of nanoblisters in the graphene sheet. These 
structural features are the origin of a very high defect density present 
in graphene after heating. Kelvin probe measurements reveal, that 
while pristine graphene flakes on mica show p-type doping beha- 
viour, defective graphene on dehydrated mica is heavily n-type 
doped. This transition is observed for bilayer graphene on mica as 
well. The exposure of defective graphene to ambient conditions 
causes again a transition from n-type back to p-type characteristic 
for pristine graphene on mica. By a detailed analysis of our data we 
could unravel the interplay of different mechanisms contributing to 
the charge carrier transfer to and from graphene. 

We would like to conclude by pointing out that our findings are 
not limited to the case of graphene on mica but have more general 
implications. From the data presented above it is evident that Kelvin 
probe force microscopy, Raman spectroscopy and even transport 
measurements performed in ambient conditions are not suitable to 
probe the intrinsic properties of graphene. These can only be 
accessed under the well-defined conditions in ultra high vaccum. 

Methods 

Sample preparation. Graphene was exfoliated a HOPG crystal (Momentive 
Performance Materials - USA) in ambient (room humidity typically 40%) onto a 
freshly cleaved muscovite mica substrate (MaTeck - Germany). Single layer graphene 
flakes were located using optical microscopy and unambiguously identified with 
Raman spectroscopy using a Horiba Labram {/i = 633 nm, P < 1 mW). 



Ambient characterization. Atomic force microscopy measurements in ambient have 
been performed using a Veeco Dimension 3100 system and Nanosensors NCHR 
cantilevers (fo ~ 300 kHz). For the Raman mappings a Renishaw InVia Raman 
spectrometer has been used. The step width has been set to 100 nm, the spotsize of the 
laserbeam is below 1 fim (/I = 532 nm, P = 0.4 mW). 

UHV characterization. In situ NC-AFM measurements (p^, = 1 X 10"^° mbar) were 
conducted using a RHK UHV 7500 system and the PLL Pro 2 controller. Highly 
conducting Si cantilevers (Vistaprobe T300) with typical resonant frequencies of/~ 
300 kHz are used for in situ AFM. Kelvin probe force microscopy^^ was performed by 
applying an AC voltage of U^c~ 1 Vand/^c^ 1 kHz. The built in lock- in amplifier 
of the PLL Pro2 controller is used to apply a DC voltage in order to minimize the 
resulting electrostatic force caused by the contact potential difference between the tip 
and surface. During the heating the temperature of the samples was monitored using a 
thermocouple. 
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